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Abstract 

We present a new approach for constructing optical phase-space-time-frequency tomography 
(OPSTFT) of an optical wave held. This tomography can be measured by using a novel four- 
window optical imaging system based on two local oscillator fields balanced heterodyne detection. 
The OPSTFT is a Wigner distribution function of two independent Fourier Transform pairs, i.e., 
phase-space and time-frequency. From its theoretical and experimental aspects, it can provide 
information of position, momentum, time and frequency of a spatial light field with precision 
beyond the uncertainty principle. Besides the distributions of x — p and t — u, the OPSTFT 
can provide four other distributions such asx — t, p — t, x — u; and p — oj. We will simulate the 
OPSTFT for a light field obscured by a wire and a single-line absorption filter. We believe that 
the four-window system can provide spatial and temporal properties of a wave field for quantum 
image processing and biophotonics. 
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INTRODUCTION 



Most of optical imaging methods are limited by the temporal and spatial resolutions be- 
cause of the unwanted scattering light coupled into their detection systems. This limitation 
is due to the fundamental concept in the process of measurement, that is, the uncertainty 
principle. Better resolution in position will reduce the resolution in momentum (angle) of an 
optical imaging system. Similarly, better resolution in time domain will reduce the flexibil- 
ities of spectroscopic analysis on a physical object. According to the uncertainty principle, 
the position x and momentum p, time t and frequency of a spatial light field cannot be 
measured simultaneously with high resolution. However, the distribution of x and p, t and 
Lo of the spatial light field can be measured simultaneously with high resolution by using 
two local oscillator fields. The use of two local oscillator fields in a balanced heterodyne 
detection scheme is also called two- window technique [H [2]. The optical phase-space to- 
mography (OPST) [H |3] is Wigner distribution [51, l5] associated with a Fourier Transform 
pair of position and momentum (angle) coordinates. Wigner distribution associated with 
two independent Fourier Transform pairs such as position-momentum and time-frequency is 
called optical phase-space-time-frequency tomography (OPSTFT) or Wigner phase-space- 
time-frequency distribution. In this paper, we develop a four-window heterodyne detection 
scheme based on two local oscillator (LO) fields for measuring the OPSTFT. The OPSTFT, 
W(x,p, a;,t), will offer the correlation information of x, p, t, and w of a wave field through 
the distributions of a; — p, u — t, x — t, p — t, x — u, and p ~ uj, where the two variables are 
plotted by fixing the other two variables. 

In quantum optics, Wigner function is usually used to represent the quantum mechanical 
wave function or quantum state of a physical system because there is no quantum device 
can directly measure the wave function. Raymer [614T0] has pioneered the measurement of 
Wigner distributions for quadrature-field amplitude of non-classical state by using optical 
homodyne detection. The method involves tomographic inversion (Radon transform) of a set 
of measured probability distributions of quadrature amplitudes. The measurement method 
developed by Raymer has also been used to measure the Wigner distribution for transverse 
spatial state (spatial mode) in a single photon level [10] and time-frequency domain [11] of 
an optical electromagnetic field. 

Recently, residual spatial fluctuation in terms of small displacements and tilts (momen- 



2 



turn) of a whole optical beam has been observed and used to exhibit EPR entanglement [12] . 
Quantum imaging [T3] has played a central role in understanding spatial fluctuations in 
quantum regions. Controlling these quantum fluctuations can improve image resolution and 
beam positioning for targeting technology. Entanglement with a large number of modes 
such as images has been accomplished by using four- wave mixing in an atomic vapor [Ti] . 
The property of multimode squeezed light [15] which allows us to increase the sensitivity 
beyond the shot-noise limit, could create many interesting new applications in optical imag- 
ing, high-precision optical measurements, optical communications and optical information 
processing. The Wigner function contains sub-Planck phase-space structures [TBI HI] ) which 
can be used to detect small movement of a large object. 

The main advantage of applying Wigner distribution in optical imaging and sensing is cer- 
tain wave-particle features of Wigner distributions, related to optical coherence, can survive 
over distances that are large compared to the transport mean free path [18]. Another ad- 
vantage is that Wigner phase-space distributions could bridge the gap between phenomeno- 
logical transport equations and rigorous wave equation treatments. Since rigorous transport 
equations can be derived for Wigner distributions, they are essential for obtaining fundamen- 
tal new insights into the nature of light propagation in multiple scattering media. Evolution 
equations for Wigner distributions, which include optical coherence scatterings, are gener- 
ally non-local and are relatively unexplored. With suitable approximations, these non-local 
equations reduce to the usual radiative transport equations [19J. Establishing the physical 
relationship between Wigner distributions and the phenomenological specific intensity will 
impact most existing methods of imaging in multiple scattering media. 

For a wave field varying in one spatial dimension and one spectra domain, S{x,u), the 
Wigner phase-space-time-frequency distribution is given by, 

W(x,p,c.,t) = / ^e- / |^e^™(r(x+ + |,c.- ^)) (1) 

where (...) denotes a statistical average. It is easy to show that J dp J dtW{x,p,ijj,t) = 
\£{x,uj)\'^, J dx J dujyV{x,p,uj,t) = \E{p,t)\'^, and etc. Most important, Eq. [l] shows that 
the Wigner distribution is Fourier Transform related to the two-point mutual coherence 
function in position and spectra. Therefore, it is sensitive to the spatially and spectrally 
varying phase and amplitude of the field. 
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FIG. 1: (a) the x-p and (b)u;-t representations of the LO fields. The shaped area is position- 
momentum and time-frequency resolutions showing beyond the uncertainty limits. 

THEORETICAL APPROACH: MEASUREMENT OF OPSTFT 

In this work, we develop a four-window heterodyne detection scheme based on two local 
oscillator (LO) fields for measuring the OPSTFT. As shown in Fig. 1(a) and (b), a local 
oscillator (LO) field is a phase coherent superposition of two fields i.e a focused and a 
collimated Gaussian fields. The focused LO Gaussian beam has spatial width of Axa = 
a and broad optical spectrum with bandwidth of Aua = a. The collimated (large) LO 
Gaussian beam has spatial width of Ax a = ^ and narrow optical spectrum with bandwidth 
of Aua = P- The position resolution a is provided by the focused LO Gaussian beam. 
The momentum resolution 1/A is provided by the collimated LO Gaussian beam. The 
purpose of this arrangement is to obtain independent control of position and momentum 
(angle) resolution such that the product of {Axa ■ Apa) = a ■ 1/A < 1 clearly surpasses 
the uncertainty principle limit as shown in the shaped area in Fig. 1(a). Simultaneously, 
the method can provide independent control of frequency (spectra) and time (path-length) 
resolution such that {Aua ■ Ata) = (3 ■ l/a < 1 to surpass the uncertainty principle limit as 
shown in the shaped area in Fig. 1(b). The spectra-resolved resolution P is provided by the 
narrowband collimated LO Gaussian beam. The path-resolved resolution 1/a is provided 
by the broadband focused LO Gaussian beam. In other words, we use position (time) 
window to obtain position (temporal) information of an optical field and simultaneously use 
momentum (frequency) window to reject noises due to other propagating angles (frequencies) 
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FIG. 2: Four-window balanced heterodyne detection for measuring OPSTFD. 

from coupling into the detection system. 

A schematic setup of the four-window technique is shown in Fig. 2. In this scheme, opti- 
cal phase-space-time-frequency tomography (OPSTFT) is measured by scanning transverse 
position (dx) and transverse momentum (dp), and by tuning wavelength (rf^) and path delay 
(r). Our system employs balanced heterodyne detection of the probe signal field, which is 
overlapped with two strong local oscillator fields (LO). The beat amplitude Vb is determined 
by the spatial and spectra overlapping of two local oscillator (LO) and signal (S) fields at 
the plane of the detector {Z = Zd) as, 

Vb = J dx' J dujElo{x', uj, zd)Es{x', u, zd) (2) 

where Elo{S){x' ,u, zd) is used to represent the two local oscillator fields (signal field), re- 
spectively. The x' denotes the transverse position at the detector plane. As shown in Fig. 2, 
when the LO fields are translated off-axis for a distance of dx and the tunable filter is 
tuned d^ away from optical center frequency, the LO field will have its spatial and spectra 
arguments shifted in Eq. |2] as given by, 

Vsidx, dui) = j dx j duElo{x' - dx,uj - d^, zd)Es{x', u, zd). (3) 

Note that the center frequency of the collimated LO beam is needed to be tuned because we 
assume the focussed LO beam has a broad optical spectrum. Using the Fresnel approxima- 
tion and standard Fourier optics technique, we can relate the fields at the detector plane, 
E{x',u, Zd), to the fields at the source planes, E{x,u, z = 0), of lenses LI and L2 as follow; 
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(i) The LO and signal fields each will experience a spatially varying phase of e 2/ after 
passing through the lens, (ii) Since the lens L2 in the signal beam is translated off-axis by 
dp and the path-length of the signal beam is delayed by r, the signal field will experience 
additional phase-shift of e 27 and delay of e~*^^, respectively, (iii) The LO and signal 
fields propagate to the detector plane through a distance of d = f, which is the focus- length 
of lenses LI and L2, and will experience the phase-shift of e 27 From (i), (ii) and 
(iii), the LO and signal fields in Eq. [sjcan be rewritten in terms of the input field z = as 
described in [2], 



Eioix' - dx,u - d^^zo) = J J dxe 2/ e ' Eioix - d^^uj - d^, z = 0), (4) 



Es{x',u,zn) = J^J dxe'^^e-'^^e-'^-Es{x,u,z = 0). (5) 
By using simple algebra, one obtains the mean square beat amplitude as given by, 

//■ . kdpX 

du / dxElo^x - dx,u - d^)Es{x,u)e~'~e~''^'' 

/r .kdpx' . , 

duo' / dx'ELo{x - 4, oj' - d^)E*s{x\ uj')e^e''^ " (6) 

Then, using the following variable transformations, a: = Xo + f,w = cJo + ^,a;' = Xo — f, and 
Lo' = Uo — \ 1 where the Jacobian of this transformation is 1. The mean square heterodyne 
beat signal of Eq. (6) can be rewritten in terms of these variables as, 
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\Vb\^ = j duodr]^dxodr]ElQ{xo - (i^; + |, Wo - 4 + ^)Elo{xo - 4 - |, - 4 - 

X E*(xo - |c^o - '^)Es{xo + |,u;o + '^)e-'"-T e-^^-\ (7) 

Recall that the Wigner function is the Fourier Transform of the two-point coherence function. 
Thus, the inverse transform for the signal field is given by, 

E*s{xo + |,u;o + '^)Es{xo " " y) = / dpdte-''^^e-'^-'Ws{xo,P,ujo,t). (8) 

By substituting Eq.(8) into Eq.(7), the mean square heterodyne beat signal is then rewritten 
by, 



duJodxodpdtdr]^dr]ElQ{xo - d^ + -,uJo - d^ + -^)Elo{xo - 4 - 77, - 4 - 



X >V(xo,p,a;o,t)e~"'(^+^)e-*^"("+*\ (9) 
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Using a similar procedure for the Wigner function of the LO field as given by, 
>Vlo(xo - d,,p + k^,00o -d^,t-T) = J dr^dr^J'^^^^'^h'^'-it - r) 

X ElQ{xo-d.^,p + k^,ujo-d^ + ^) 

X Elo{xo - d^ - - - Uo - d^ - y). (10) 

Then, the mean square heterodyne beat signal is finally obtained as given by, 

|VbP = Jdxjdu JdpJ dtWLoix - d^,p + k^,uj - d^,t - T)Wsix,p,uj,t)) (11) 

where we have changed the notations of Xo ^ x and Uo — )■ uj. The |VbP is proportional to 
the phase-space-frequency-time convolution integral of the Wigner distributions for the two 
local oscillator and the signal fields in the planes of the input lenses LI and L2, respectively. 

We will use a broadband light source for this experiment. The LO fields are engineered 
in the form as given by, 

EUx^oo) = ^o[exp(-^)exp(-^)) +7exp(-^)exp(-^)e^^] (12) 

where the spectrum of LO fields are assumed to be Gaussian function. This can be accom- 
plished by using a single mode fiber with a tunable bandpass Gaussian filter from Newport. 
The phase- dependent part of Wigner function for the LO takes the form 

2x^ 

WLo{x,p,uj,t) (X exp[— -2aVH ^ - 2/3^^^] cos(2xp + 2wt + 0) 

^ cos{2xp + 2ujt + (p) (13) 

where we take the approximation of A ^ a and a ^ (3. Then, the range of integration 
for the momentum, position, frequency and time coordinates in Eq. [Tl] is limited by the 
signal field. In this scheme (similar electronic components as in ref[ [HE]]), the two LO 
fields which differ in frequency by kHz are phase-locked. The signal field is modulated 
such that the heterodyne beat signals with the focused LO field and the collimated LO 
field are about MHz and MHz + kHz, respectively. The root mean square beat 
amplitude is measured with an analog spectrum analyzer with a bandwidth of 100 kHz (> 
kHz) centered at fi^ MHz. The output of the spectrum analyzer is squared in real time 
with a low noise amplifier and an analog multiplier, then the amplified signal is sent to the 



lock-in-amplifier. Substituting Eq. 13 into Eq. 11, we find that the in- and out-of phase 



quadrature amplitudes in the lock-in-amplifier are directly corresponding to the real and 
imaginary parts of the quantity, 



|VbP oc }C{Xo,Po,U}o,to) 

oc {Eg{ )Es{po, to))exp{iXoPo + iuJoto) 

oc SR + iSi (14) 



The /C(x, p, CO, t) is the Kirkwood-Rihaczek phase-space-time-frequency distribution. Eq. 14 
is readily inverted to yield the Wigner phase-space and time-frequency function or the OP- 
STFT of the signal field by a linear transformation. We obtain, 

Ws{x,p,uj,t) oc J dxodpodujodto cos[2(a: - Xo){p - Po) + 2(a; - 0Jo){t - to)]SR{xo,Po,(^o,to) 

+ J dXodpoduJodtoSm[2{x - Xo)ip - Po) + 2{uj - uJo)it - to)]Sj{xo,Po,uJo,to) 

(15) 

which is the OPSTFT of the signal field. Sr and Sj are the real and imaginary parts 



of Eq. 14, i.e. the in- and out-of-phase quadrature amplitudes, which are simultaneously 



measured. 



SIMULATION OF OPSTFT 
A Thin Wire 

The measurement of OPSTFT based on two local oscillator fields is called four-window 
technique. In practice, the four independent variables a, I /A, 1/a and (3 are chosen to be 
small to resolve scales of interest on the spatial and temporal properties of a wave field. 
Physical properties of an object can be extracted through measuring the OPSTFT of the 
scattered light field through the object. The OPSTFT can provide x — p and u — t distribu- 
tions, and other four distributions such a.s u — x, t — x, u — p and t —p distributions. These 
six distributions can provide new types of information for the light field under study. We 
numerically simulate the measurement of OPSTFT for the light field scattered through a 
thin wire with the diameter of 0.6 mm. We use a Gaussian beam with a wave field as given 
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FIG. 3: The Kirkwood-Rihaczek (/C) distribution (a) real and (b) imaginary parts of /C(x,p, 0, 0). 
(c) real and (d)imaginary parts of /C(0.4, 0, c^;, t) 



by, 



X 



S{x,t) oc exp[-^]exp[2^]exp[-^]exp[za;oto], (16) 

where and at are the spatial and temporal bandwidths. R and Uo are the radius of 
curvature and center/carrier frequency of the light field. We are interested in looking at the 
light scattered right after the wire, where the scattered light field, Swire{x, t), can be written 



as the product of Eq. 16 and a wire function (slitfun[x]=If[— 0.3mm < a; < 0.3mm, 0.0, 1.0] 
in Mathematica program). We have used this field for exploring phase-space interference 
analog to superposition of two spatially separated coherent states [2]. However, the phase- 
space distribution associated with time-frequency distribution incorporated into Wigner 
distributions haven't been explored. We use cr^. = 0.85 mm, at=200 fs, and -R=-10000 
mm and Uo = for the simulation. We first obtain Smre{p,^) by numerically Fourier 
transforming the Sujire{x,t). Then, we generate the Kirkwood-Richaczek phase-space-time- 



frequency distribution, )C{x,p,u,t) = t)£^i„ire(p, i^)exp(za;p — zwt), as in Eq. 14 for the 



scattered field S^ireix, t). We first plot the real and imaginary parts of IC{x,p, 0, 0) as shown 
in Fig.3(a) and (b). We adopt the units for position (x) in mm, momentum (p) in rad/mm, 
frequency (u) in 10^^ Hz, and time in 10^^'^s. The }C{0,0, uj,t) is zero because there is zero 
beat signal (|VbP) at the position x=0 and p =0. To explore the time-frequency distribution 
associated with nonzero beat signal (|VbP) at the phase-space point(x=0.4, p=0), we plot 
the real and imaginary parts of /C(0.4, 0, w, t) as shown in Fig.3(c) and (d). From the 
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FIG. 4: The optical phase-space-time-frequency tomography (OPSTFT) (a)W(x,p, 0, 0), 
(b)>V(0,0,a;,t), (c)W(O, 2, w, i), (d)W(0,p, 0, t), (e)W(0,p, w, 0), (f) W(x,0,0,t), and 0, a;, 0) 



]C{x,p,u,t) distribution, we can obtain Wigner phase-space-time-frequency distribution by 



using linear transformation as in Eq. 15 We use Mathematica program for performing the 
numerical transformation. Fig.4(a) is the plot of VV(x,p, 0,0) where the tunable filter in 
the collimated LO beam is set at center frequency and the delay (r) in the signal beam 
is set to zero. The phase-space oscillation along the x=0 is due to the coherent phase- 
space interference of two spatially separated wave packets after the wire. Now, we plot the 
>V(0,0,a;,t) as shown in Fig.4(b), where the phase-space point is at (0, 0). Note that the 
/C(0, 0, u, t) is zero everywhere but not the Wigner function of W(0, 0, u, t). The reason is /C 
distribution is suitable for describing local properties of a wave function, while the Wigner 
function is suitable for describing wave properties of a particle wave function [2]. Fig. 4(c) 
shows an interesting result of W(0, 2, cj, t), which is negative, i.e, the inverse of Fig. 4(b). 
This is because the Wigner distribution of W(0, 2, 0, 0) has negative value. These properties 
are important to explore hyper-entanglement of intrinsic properties of single-photon spatial 
qubit states. Other four distributions such as >V(0,p, 0,t), W(0,p, a;,0), >V(x,0,0,t), and 
W{x, 0, u, 0) are plotted as shown in Fig.4(d), (e), (f) and (g), respectively. The >V(0,p, 0, t)) 
and W(0,p, u, 0) exhibit oscillation/interference behavior along momentum coordinate. The 
phase-space interference due to the two spatially separated wave packets influenced the 
distribution of momentum/angle in the time and spectra domains. 
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FIG. 5: The Kirkwood-Rihaczek (/C) distribution (a) real and (b) imaginary parts of /C(0, 0, w, t). 
(c) real and (d)imaginary parts of /C(x,p, 0.2, 0) 

An Absorption Filter 



We numerically simulate the measurement of OPSTFT for a broadband light field passing 
through an absorption filter with bandwidth of 2 THz. We use a Gaussian beam with a 
Gaussian linewidth as given by, 



S{x, io) oc exp[- 



l^jexp^.— jexp^ 



(17) 



where a^^ is the spectra bandwidth. It is much easier to work on the spectra domain of the 
field so that the light field passing through the filter can be written as the product of Eq. 17 
and an absorption filter function (slitfun[x]=If[— 0.1 < x < 0.1, 0.0, 1.0] in Mathematica 
program). We use ax = 0.85 mm, crt=5.0 THz, i?=-10000 mm and Uo = for the simula- 
tion. First, we generate the Kirkwood-Richaczek phase-space-time-frequency distribution, 
)C{x,p,u,t) = £J^i^^j.{x,uj)£fiiteT{p^'t)^^'9{.'''^'P + i<^t), where the £fiiter{p,t) is obtained by nu- 
merically Fourier transformed the Sfuterix, uj). The ]C{x,p, 0, 0) is zero because there is zero 
beat signal (|VbP) at the position u = and t = 0. Fig.5(a) and (b)show the real and imag- 
inary parts of /C(0, 0, w, t). To explore the position-momentum distribution associated with 
nonzero beat signal (|VbP) at the time-frequency point(w=0.2, t=0), we plot the real and 
imaginary parts of /C(0.2, 0, u, t) as shown in Fig. 5(c) and (d). Since we have the /C(x,p, u, t) 
distribution, then we can obtain Wigner phase-space-time-frequency distribution by using 



linear transformation as in Eq. 15 
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FIG. 6: The optical phase-space-time-frequency tomography (OPSTFT) (a)W(0, 0, w, t), 
(b)W(x,p,0,0), (c)W(x,p,0,3), (d)>V(0,p,0,t), (e)W(0,p,a;,0), (f) W{x,0,0,t), and 0, tj, 0) 

Fig. 6(a) is the plot of W(0,0,a;,t) where the position dx of the mirror for the LO beam 
is set to zero and the momentum dp of the lens LI in the signal beam is set to zero. The 
time-frequency oscillation along the a; = is due to the coherent time-frequency interfer- 
ence of two spectrally separated wave packets after the filter. This observation has been 
observed in Wigner time-frequency distribution [TT]. We plot the W(x,p, 0,0) as shown in 
Fig.6(b), where the time-frequency point is at (0, 0). Fig.6(c) shows an interesting result 
of W{x,p, 0,3), which is negative, i.e, the inverse of Fig.6(b). This is because the Wigner 
distribution of W(0, 0, 0, 3) has negative value. Other four distributions such as W(0, p, 0, t), 
W{0,p,ijj,0), W{x,0,0,t), and W(x,0,aj,0) are plotted as shown in Fig.6(d), (e), (f) and 
(g), respectively. The yV{0,p, 0, if:)) and W(x, 0, 0, t) exhibit oscillation/interference behavior 
along p =0 and x=0, respectively. 

CONCLUSION 

In conclusion, we have developed four-window optical heterodyne imaging technique 
based on two local-oscillator (LO) beams for measuring OPSTFT of a wave field. The 
four-window technique can simultaneously provide high resolution in position, momentum 
(angle), time, and spectra for characterizing spatial properties of a wave function. This 
newly developed optical phase-space-time-frequency tomography can be used for explor- 
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ing hyper-entanglement spatial qubit states of photon wave function and for early disease 
detection in biophotonics. 
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